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The effects induced by long temporal correlations of the velocity gradients on the dynam- 
ics of a flexible polymer are investigated by means of theoretical and numerical analysis 
of the Hookean and FENE dumbbell models in a random renewing flow. For Hookean 
dumbbells, we show that long temporal correlations strongly suppress the Weissenberg- 
number dependence of the power-law tail characterising the probability density function 
(PDF) of the elongation. For the FENE model, the PDF becomes bimodal, and the 
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coil stretch transition occurs through the simultaneous drop and rise of the two peaks 
associated with the coiled and stretched configurations, respectively. 



1. Introduction 

The dynamics of a flexible polymer in a moving fluid depends strongly on the properties 
of the velocity gradients. A remarkable e xample is given by the coil-stretch transition i n 
a planar extensional flow v — (ex, —ey) ( de Gennea 1974 IPerkins. Smith fc Chu 1997 ). 
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I This flow is characterised by a direction of pure compression and a direction of pure 

Q\ . stretching with constant strain rate e. The elongation of the polymer is controlled by the 

Weissenberg number Wi = er p , where t p is the longest relaxation time of the polymer. 
. The value Wi = 1/2 marks the coil-stretch transition: for Wi < 1/2 the polymer remains 

coiled under the action of the entropic elastic force, whereas for Wi > 1/2 the drag 
force exerted by the flow overcomes the entropic force and the polymer unravels almost 
completely. As a consequence, the probability density of the extension, p(R), has a single 
pronounced peak, whose position depends on Wi. For Wi < 1/2, the peak is in the 
neighbourhood of the equilibrium size of the polymer, i?o; for Wi > 1/2, the peak 
approaches the maximum length i? m ax- 

The ability of a (non-uniform) laminar flow to deform an isolated poly mer has now 
bee n demonstrated for various flow configurations (see, e.g., the reviews by lLarson 120051 



and ISha qfeh l2005h. The corresponding problem for random flows was first studied by 
iLumlevi (jl97aT Il973l ). who observed that the strain tensor of an incompressible 



ran- 
dom flow always has a direction of stretching, although such direction fluctuates in time 
and in space. Lumley further remarked that the vorticity penalises polymer stretch- 
ing, for it prevents the polymer from remaining aligned with the principal axes of the 
strain tensor. He then concluded that random flows can stretch polymers if the product 
of the Lagrangian correlation time of the strain tensor and the modulus of the maxi- 
mum eigenvalue of the deformation tensor exceeds a critical value proportional to r" 1 . 
iGroisman &; Steinberg! (|200ll ) confirmed these conclusions by showing evidences of a sig- 
nificant amount of polymer stretching in a low-Reynolds-number random flow generated 
by viscoelastic instabilities. 

In random flows, the time scale describing the stretching of line elements is the re- 
ciprocal of the maximum Lyapuno v exponent A. Hence the appropriate d efinition o f the 
Weissenberg number is Wi = Ar p . iBalkovskv. Fouxon fc Lebedev ( 2000h and Chertkov 



2 S. Musacchio and D. Vincenzi 

(2000) related the deformation of a polymer in a random flow to the statistics of the 
stretching rate or, more pre cisely, to the entropy function ass ociated with its proba- 
bili ty density function (e. g., Crisanti. Paladin fe Vulpiani 19931) . Based on this analy- 
sis, iBalkovskv et al. I (<2000l) predicted the existence of the coil-stretch transition for any 
random flow with positive A, and proposed the following explanation. For an infinitely 
extensible polymer (i? max — > oo), the stationary probability density function (PDF) of 
the extension has a power-law tail: p(R) oc i? _1 ~ Q for R ^S> Rq, where a depends on 
the form of the entropy function. The exponent a is positive for Wi < 1/2, it decreases 
with increasing Wi, and becomes negative as Wi exceeds 1/2. Consequently, p(R) is 
normalizable only for Wi < 1/2, whereas for Wi 1/2 there is unbounded growth of the 
extension and the assumption of infinite extensibility becomes inadequate. This abrupt 
change in the statistics of the extension is interpreted as indicating the coil-stretch tran- 
sition in random flows. The threshold Wi — 1/2 reproduces (to within a numerical factor) 
Lumley's criteriorjj 

If the random flow is isotropic and has zero correlation time, p(R) can be written 
explicitly and the behaviour predicted by Balkovskv et al.l ( 2000l) can be obtained by 
direct computation ( Chcrtko vl2000t lThiffeaultll2003t ICelani. Musacchio fc Vincenzill2005t 



Martins Afonso fc Vincenzil 120051 ) . In this particular case, a = d(Wi — 2)/2, where d 



is the spatial dimension of the flow. 

A polymer having finite maximum extension (i? m ax < oo) reaches a stationary configu- 
ration at any Wi since the entropic force forbids extensions greater than i? m ax- According 
to the analysis of the infinitely extensible case, p(R) is now expected to display a power- 
law behaviour R~ 1 ~~ a for Rq <C R <C i?max- In a short-correlated flow, the coil-stretch 
transition results from the combination of this power law and the cutoff at R = R ma , x : 
as Wi increases, the interm ediate power law raises and the pe ak of p{R) moves from Rq 
to extensions near to i? max ([Martins Afonso fc Vincenzill2005l ). An important qualitative 
difference between extensional and random flows must nevertheless be emphasised. In 
the former case, p(R) has a pronounced peak at an extension near either Rq or i? max 
depending on whether Wi is less or greater than 1/2. In the latter case, p(R) has broad 
tails signaling the coexistence of coiled and stretched configurations with predominance 
of either configuration according to the value of Wi. In random flows, the coil-stretch 
transition is not as sharp as in the extensional case because of the random nature of the 
velocity gradient and the presence of vorticity. 

The coil-stretch transition and the relation between th e exponent a and the maximum 

Lyapu nov e xponent A were investig ated experimentally bv lGerashchenko. Chevallard fc Steinberg 
( 20051 ) and Liu fc Steinberg! ( 2010t ). In these experiments, the statistics of the extension 
was measured by following a fluorescently labelled polymer in a random velocity field 
generated by elastic turbulence. The PDF of polymer extension was also investigated by 
means of direct numerical simulations of the continuum equation for the polymer confor- 
matio n tensor ( Eckhardt. Kroniager fc Schumacherl l2002: Bof fetta. Celani fc Musacchio! 

20031) and by means of Brownian Dynamics simulations for elastic dumbbells ( Celani et al.l 

2005 [lDavoudi fc Schumacherl200a[Bagheri et al.l2010h or multi-bead chains (Watanabc fc Gotohl 
20101) . The results of the simulations support the picture of polymer dynamics presented 



above. 

The Weissenberg number suffices to determine the elongation of a polymer in lami- 
nar flows. By contrast, in random flows an additional dimensionless number may influ- 
ence the deformation of the polymer, namely, the Kubo number Ku = Ar c , where r c 



f IBalkovskv et al.l ()2000l ) defined the Weissenberg number as Wi = 2\t p so that the 
coil-stretch transition occurs at Wi = 1. 
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is the Lagrangian correlation time of the velocity gradient. The case Ku — has been 
briefly reviewed above. In this paper, we examine the effect of a nonzero Ku and show 
that the dynamics at Ku > 1 differs significantly from the one predicted for short- 
correlated flows. To emphasise the basic physical mechanisms, we co nsider a simplifie d 



situation. The polymer molecule is modelled as an elastic dumbbell (jBird et al.l 119771 ). 



This app roximation is appropriate when attention is restricted to the extension of the 
polymer ( Watanabe fc Gotoh 201Clh . As for the random advection, we consider a two- 
dimensional linear renewing flow. In renewing flows, time is split into intervals of length r c 
and the velocities in differ ent intervals are independent and identically distributed (e.g., 
Child ress fc Gilbertlll995l p. 320). Linear renewing flows were used by IZeT dovich et al.l 
(| 1984T ) to study the kinematic dynamo effect. The magnetic field is actually stretched 
by the velocity gradient in the same way as an infinitely extensible polymer. The use 
of a linear renewing flow enables us to obtain semi-analytical results and to accurately 
compute the statistics of polymer extension with moderate numerical effort. 

The remainder of the paper is divided as follows. In section [2j we briefly review the 
elastic dumbbell model. In section[3J we introduce the renewing random flow. The results 
are presented in section [4] Finally, some conclusions are drawn in section [5] 



2. Elastic dumbbell model 

An elastic dumbbell is composed of two beads joined by a spring. The beads repre- 
sent the two ends of the polymer; the spring models the entropic force. For the sake of 
simplicity, we assume that the flow transporting the dumbbell is two-dimensional. 

The vector separating the ends of the polymer, R, satisfies the stochastic ordinary 
differential equation (e.g., Bird et al. 19771 ): 



R = a (t)R-^-f(R)R+J^m (2.1) 

where R = \R\ and £(i) is white noise, i.e., a Gaussian process with zero mean and 
two-time correlation (t;i(t)£,j(t')} — 5ij8(t — t'). The three terms on the right-hand-side 
of equation (|2.1[) result from the drag force, the entropic elastic force, and thermal noise, 
respectively. The 2x2 matrix a(t) is the velocity gradient evaluated along the trajectory of 
the centre of mass of the dumbbell: <Ty (i) = djVi(t). The function f(R) is identically equal 
to 1 for an infinitely extensible polymer (Hookean model) and has the form f(R) — 1/(1 — 
R 2 /R 2 nllx ) for a finitely extensible polymer with nonlinear elasticity (FENE model). In 
the former case, R is defined on M. 2 ; in the latter case, R belongs to [0, i? max ) x [0, -Rmax)- 
Equation (I2.1|) holds under some assumptions on the dynamics of the beads. The 
velocity field is assumed to be linear at the size of the dumbbell. The drag on a bead 
is given by the Stokes law. Furthermore, hydrodynamic interactions between the beads 
and inertial effects are disregarded. 



3. Linear renewing flow 

According to the assumptions of the dumbbell model, in the reference frame of the 
centre of mass of the polymer the velocity field is of the form: Vi(r,t) = aij(t)rj. The 
velocity gradient o~(t) is constant in each of the time intervals /„ = [nr c , (n+ l)r c ), n G N. 
We denote by er„ the (random) value taken by a(t) in /„: a{t) = a n for all t £ I n . The 
random matrices a n are identically distributed and statistically independent. We assume 
that, for a fixed n, a n is Gaussian, zero-mean, statistically isotropic, and traceless (so 
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that the velocity field is incompressible). As a result, a(t) takes the form: 

with S and J? positive constants. The elements CtW satisfy: Ci(^) = Ci,n f° r an i 6 4i 
where the Ci,n are Gaussian random variables such that (Ci,n) = for all i and n and 
((i,n(j,m) = SijS nm for all i, j — 1,2,3 and n, m G N. The mean-squared strain and 
rotation rates are S 2 and .f? 2 , respectively. We set J? = 5* to reprodu ce the relation 
holding for the solution of the Navier-Stokes equation (e.g.. iFrisch 19951 p. 20). 



Strictly speaking, the flow considered is not statistically stationary in time. Neverthe- 
less, over time ranges longer than r c , it can be considered as a good approximation to 
a stationary flow with correlation time t c owing to the invariance with respect to the 
transformation t — » t + r r and thanks to the independence of the matrices a(t) in different 



time intervals /„ (jZerdovich et al.lll984l) . The <5-correlated flow (i.e. white-in-time noise) 
is recovered by letting r c tend to zero while holding S 2 t c constant. 

As mentioned in the introduction, the elongation of a polymer is related to the statis- 
tics of the stretching rate of the flow. For a review on the entropy function and the 
ge neralized Lyapunov ex ponents in statistical physics, we refer the reader to the book 
bv ICrisanti et ah (1993). Here, we briefly recall some basic elements of the theory. If 



l(t) denotes a fluid-line element, the stretching rate at time t is defined as 7 (i) 
t~ x hi[£(t)/£(0)}. The maximum Lyapunov exponent is the long-time limit of the av- 
erage stretching rate: A = lim t _ ) . 00 (7(i)), where the average is ta ken over the realizati ons 
of a(t). The PDF of 7 for long t takes the large-deviation form ( Crisanti et al. 1993h : 



P( 7 ,t)oce- G(7)t , 

where G( 7 ) is the entropy function. G{^) is non-negative and attains its minimum value 
at the point A. Equivalently, the stretching properties of the flow can be characterised 
by the generalized Lyapunov exponents, defined as the rate of growth of the moments 
of£(t): 

L(,)=limiln/fMy 
W t^oot \\t(Q)J 

The function L(q) is the Legendre transform of G( 7 ), L(q) = max-j^g — 6(7)], and sat- 
isfies -Zj'(O) = A. For small q, L(q) can be obtained by using the quadratic approximation 
G( 7 ) « (7 - A) 2 /(2A) so that L(q) « Xq + Aq 2 /2. 

For the linear renewing flow, we can compute L( q) by adap ting to the present problem 
the method described by iGilbert fc BavTvl (19921) (see also IChildress fc Gilbert! Il995l 



pp. 322-326). A line element satisfies the equation: 

£(t) = a(t)£(t). (3.2) 

Since a(t) is constant in I n , the solution of equation (|3.2j) for t € I n is: l(t) = e cr " t £(nr c ). 
Hence 

£((n + l)r c ) = J n £(nr c ), 

where J n = e anTc is written 

2w n c+ + SCi {St 2 ,n + V2f2( 3 



1 

2uj t . 



\(S( 2 ,n - V2ft(3,n)c n 2w„C+ ~ S(l in C r , 



with w n = US*(g n + a n ) - 2^ 2 C 3 2 „, 4 = (e^ + e-»^)/2, and c" = (e«^ 
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Figure 1 . Lyapunov exponent as a function of t c (circles) . Dashed and dotted lines represent the 
asymptotic behaviours (|3.5j) for r c — > and r c — > oo, respectively. Inset: Generalized Lyapunov 
exponents L(q) for three values of Ku. Note that L{— d) = L(0) = in ag r eemen t with the 
properties of the generalized Lyapunov exponents proved bv lZel'dovich et al.l (1 19841 ). 



e~ UnT ")/2 (we remind the reader that Ci,n is the constant value taken by the random 
function Q (i) in I n ). If uj n is real, then J n can be written in terms of hyperbolic functions; 
for a non-real complex lj„, J n involves trigonometric functions. By using the statistical 
isotropy of the flow and the fact that the matrices a n are identically distributed and 
indep endent for different n, it is possible to show the following relation ((Gilbert fc Bavlv 
Childress fc GilbertJll995h : 



mdep 
1992; 



(e«((n + l)T c )) 
(l"(nT c )) 



(\Jne\ q ) In 



(3.3) 



where e is any constant unit vector and the average (•)/ n is taken over the statistics of 
the gradient in the interval only. Iterating equation (13.31) yields: (l q (nT)) = (p q ) r } £ q (0) 
with p 2 = | J n e\ 2 — e T J^J n e. It follows that L(q) and A can be written as 



Hp" 



In 



and 



A = t c 1 (lnp)j 



(3.4) 



Given that the o~ n are identically distributed, the average can be taken over any in- 
terval /„. Equations (|3.4p provide a simple way to compute the function L(q) and the 
Lyapunov exponent. The behaviour of A as a function of r e is reported in figure [TJ The 
following asymptotic behaviours hold (jChertkov et al. I ll99fl) : 



A 



S 2 r c 



(t c —> 0) and A ~ Re(o; n )/ n 



7T S 

2 2 



S 2 



{ 7~c y OO j 



S 2 + 2Q 2 

(3.5) 

For Q = S, one obtains A ~ yj n/6 S/2 as t c — > oo. Thus the Lyapunov exponent 
becomes independent of r c for large r c , and the convergence to the asymptotic value is 
exponential. By contras t, the variance A mon otonically increases like A ~ S 2 t c /A both 
for small and large t c ( Chertkov et al. 199(f), and therefore the generalized Lyapunov 
exponents do not saturate to a constant value (figure Q]). We stress the fact that the above 
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Figure 2. PDFs of polymer elongation, p(R), for the Hookean model at Wi = 1/4. Note the 
power-law behaviours p(R) oc R d ~ 1 for R <C Ro and p(R) oc R" 1 ^" for i? S> -Ro. Inset: Exponent 
q as a function of .K«. Here ao = lim_ff u ^o a. Dash-dotted and dashed lines are the asymptotic 
behaviours (|4.1[1 and (|4.2I) for Ku — > and /Tit — >• oo, respectively. 



asymptotic behaviours hold for any two-dimensional random flow that is incompressible 
and statistically isotropic and not only for the renewing flow. 



4. Statistics of polymer extension 

In order to investigate the influence of temporal correlations of the velocity gradients 
on the dynamics of polymers, we numerically integrated equation (12. ip for the elastic 
dumbbell model, where the velocity gradients were determined by equation (|3.ip . We 
computed the PDF of the polymer elongation R for various values of Ku and Wi. The 
nu merical integrat ion has been performed using the predictor-corrector scheme proposed 
bv lOttingerl (1996) fo r the F ENE model and the stochastic Runge-Kutta algorithm intro- 
duced bv lHonevcuttl (|1992T ) for the Hookean model. The statistics of polymer elongation 
has been computed by following the dynamics of a single dumbbell for 10 time inter- 
vals /„. 

Within the Hookean model, t he PDF of the extension behaves like p(R) oc Rr 1 ^" for 
R > R (|Balkovskv et al J 12000). The coil-stretch transition is identified by the change 



of sign of a, which occurs at Wi = 1/2. The effects of the temporal correlation of the 
velocity gradients can be quantified through the dependence of a on the Kubo number. 
The numerical simulations of the Hookean model show that for small values of Ku the 
tail of p(R) is almost independent of Ku; conversely, for Ku > 1, the tail raises as Ku 
increases and approaches the slope —1 (figure [2]). Because of long temporal correlations, 
even polymers with a short relaxation time occasionally experience significant stretching 
events, and this effect produces a power-law tail close to the coil-stretch transition also 
for small values of Wi. 

This intuitive idea can be rationalized by the following argument. The e xponent a de- 



term ines the lowest order such that (R a ) diverges and satisfies a = 2r p L(a) (Boff etta et al 



20031 ). For Wi near to 1/2, a is not far from zero and it is appropriate to use the quadratic 
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approximation for L(q). The equation a = 2t p L(o) then yields ( Balkovskv et al. 2000h 



a = A ( J_ - 2 V (4.1) 
A V Wi J v ' 

While the dependence of a on the Weissenberg number is contained entirely in the second 
factor, the Kubo number enters the expression for a through the ratio A/A. For the linear 
renewing flow, the asymptotic behaviours discussed in section [3] give: 

A/A - 1 (Ku-> 0) and A/A = O^u" 1 ) (Ku -> oo). (4.2) 

One obtains that for small values of Ku the tail of p(R) is independent of Ku, whereas 
a = (^(Ku^ 1 ) as Ku is increased, in agreement with our numerical findings. 

Although for an infinitely extensible polymer p(R) is normalizable only for a > 
(equivalently for Wi > 1/2), the calculation leading to the power-law behaviour remains 
formally valid also if a < (Wi > 1/2). In this case, the stationary PDF of the extension 
exists only if the nonlinearity of the elastic force is taken into account, and the power-law 
prediction is expected to hold for intermediate extensions Rq -C R <C i?max- To allow the 
development of the intermediate power law, we performed numerical simulations of the 
FENE model at artificially high R max /Ro- For Wi < 1/2 the results agree with those 
obtained for the Hookean model (figure [3l panel (a)). For Wi > 1/2 an increase of Ku 
produces a decrease of the intermediate slope of p(R), which can be negative even for 
very large Wi (figureEJ panel (b)). According to equation (|4.2j) . near Wi = 1/2 the range 
of variation of a as a function of Wi can be made arbitrarily small by increasing Ku, 
to the extent that the intermediate slope of p(R) becomes almost independent of Wi 
(figure S]) . Our findings show that the picture of a coil-stretch transition associated with 
the increase of the intermediate slope of p(R) does not hold in a long-correlated flow. At 
large Ku, the appearance of the stretched state occurs through the drop of the maximum 
at R Rq and the simultaneous rise of a second maximum near i? ma x (figure |4|). 

For realistic values of Rmax/Ro, it is difficult to detect a clean power law at inter- 
mediate extensions and the above scenario becomes even more relevant. The numerical 
simulations of the FENE model for R max /Ro = 20 show that, for short-correlated flows 
(Ku <1), the coil-stretch transition manifests through the gradual widening of p(R) and 
the displacement of its maximum from values close to Rq towards values close to i? max 
(see figure [3 panel (a)). This beh aviour is in agreement with the t heoretical prediction 



for the (S-correlated isotropic flow ([Martins Afonso fc Vincenzill2005l ) an d with recent ex- 



perimental measurements of polymer elongation in a random shear flow (jLiu fc Steinberg 
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Figure 4. PDFs of polymer elongation for the FENE model at various values of Wi in a 
long-correlated flow (Ku — 10). The maximum extension is set to 7? max = 10 3 i?o- 

2010f h For long-correlated flows (Ku > 1), the dependence of p(R) on Wi is very differ- 
ent. In this case, increasing Wi at fixed Ku rather produces the drop of the peak at Rq 
and the simultaneous formation of a second peak near i? max (see figure^ panel (b)). The 
PDF of the extension is bimodal and the two maxima are clearly distinct: intermediate 
extensions are much less likely than in the small- Ku case. 

The shape of p(R) for large Ku can be explained intuitively as follows. If the Lagrangian 
correlation time of the velocity gradient is long, a typical trajectory of a polymer is 
composed of long portions where the gradient can be thought as frozen with a fixed 
direction of stretching. Some of those portions will be characterised by a weak stretching 
rate, some by a strong one. The response of the polymer to the velocity gradient along 
one of those portions of the trajectory will be similar to the response that a polymer 
would have in an extensional flow with strain rate comparable to the stretching rate. 
As mentioned earlier, the PDF of R in an extensional flow has a sharp peak either at 
small extensions or at long extensions depending on the strain rate, and intermediate 
extensions are very unlikely. Therefore, the evolution of an isolated polymer at large Ku 
consists of a sequence of long deformation events, whose intensity can produce either a 
coiled or a highly stretched extension, whence the bimodal shape of p(R) ■ 

5. Conclusions 

We have studied the influence of long temporal correlations on the dynamics of a 
Hookean dumbbell and of a FENE dumbbell in a two dimensional random renewing 
flow. The time correlation of this simple model flow can be changed arbitrarily, and this 
property enabled us to obtain semi-analytic predictions for the Kubo-number dependence 
of the statistics of polymer elongations. It is known that the PDF of the elongation of 
Hookean dumbbells is characterised by a power-law tail and that this power law also 
describes the behavior of the PDF of FENE dumbbells at intermediate elongations. We 
have shown that in the long-correlated limit (Ku —> oo) the power-law tail becomes 
almost insensitive to Wi and its slope is not an effective indicator for the coil-stretch 
transition. In the case of the FENE model, the long temporal correlations strongly affect 
the shape of the PDF, which shows two distinct maxima corresponding to the coiled and 
stretched configurations. For Ku > 1 we have found a new scenario for the coil-stretch 
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Figure 5. PDFs of polymer elongation, p(R), for the FENE model, at various Wi numbers, 
in a short-correlated flow (Ku — 1CF 2 , panel (a)) and in a long-correlated flow (Ku — 10, panel 
(6)). The maximum extension is fixed to i? max = 20i? . 



transition, which occurs through the drop of the peak at R and the simultaneous rise 
of the second peak near i? max . 

Our findings can be understood in terms of basic properties of long-correlated flows, 
and the underlying mechanisms does not rely on the peculiar choice of the model flow 
used in our study, ft is therefore arguable that the phenomena discussed here could be 
observed also in realistic flows. In particular, it would be interesting to investigate by 
means of numerical simulations or experimental measurements whether the presence of 
long-lived structures in high-Reynolds-number flows could result in the appearence of 
the two-peak coil-stretch transition scenario depicted here. 

We are grateful to Antonio Celani and Prasad Perlekar for fruitful discussions. 
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